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Abstract 

Space-time correlations between produced particles, induced by the composite nature of hadrons, imply specific changes 
in the properties of the correlation functions for identical particles. The expected magnitude of these effects is evaluated 
using the recently published blast-wave model analysis of the data for pp collisions at y/s = 7 TeV. 
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1. It has been recently pointed out [1] that since hadrons 
produced in high-energy collisions are not point-like ob¬ 
jects, they cannot be uncorrelated. Indeed, being com¬ 
posite, hadrons cannot occupy too close space-time points 
(because at small distance the constituents of hadrons mix 
and there are no separate hadrons to interfere). Conse¬ 
quently, since the HBT experiment measures the quantum 
interference between wave functions of hadrons, it cannot 
see hadrons which are too close to each other. Therefore 
the distribution function of the pair of hadrons must van¬ 
ish at small distances between them. 

This implies of course a correlation in space-time. As 
this correlation is the necessary consequence of the com¬ 
posite structure of hadrons (and thus it is a general prop¬ 
erty of the system) it is interesting to investigate to what 
extent it modifies the accepted ideas about the quantum 
interference which are, usually, derived under the assump¬ 
tion that such correlations can be neglected [2]. 

It was already shown in [1] that such space-time cor¬ 
relations may be responsible for the observation that the 
two-pion Bose-Einstein correlation function takes values 
below unity [3-5], at variance with the well-known theo¬ 
rem valid when the correlations are ignored [2]. 

In the present paper the investigation of this phe¬ 
nomenon is continued, using the recently published [6] 
analysis of the data on HBT radii, measured by the ALICE 
collaboration [7]. This allows to estimate quantitatively 
the magnitude of the effect and to give predictions for its 
size in all three directions long, side and out , commonly 
used in discussion of the quantum interference [2]. 

In the next section the consequences of the space-time 
correlations for the HBT correlation functions are ex¬ 
plained. In Sections 3 and 4 the blast-wave model used for 
the quantitative estimate of the effect is presented. The 
results are presented and summarized in the last two sec¬ 
tions. 

2. In absence of correlations between produced hadrons, 


the two-particle source function is the simple product 

w{pi,p 2 ;x 1 x 2 ) = w(p 1 ,x 1 )w(p 2 ,x 2 ) (1) 


where w{p, x ) is the single-particle source function (Wigner 
function). Consequently, the Bose-Einstein correlation 
function between the momenta of two identical particles 


C(pi,P2) 
is given by [2] 
C(pi,p 2 ) 


Here 


w(P i2 ;Q) 

w(p) 


N{pi,P2) 

N( Pl )N( P2 ) 


1 + w(Pi 2 ;Q)w(Pi 2 ;-Q) 
w(p 1 )w(p 2 ) 

t , \w(Pi2,Q)\ 2 >L 

w(pi)w(p 2 ) - 


= J dx e lQx w(P 12 ;x), 
= Jdxw(p;x), 


( 2 ) 


( 3 ) 


( 4 ) 


where P i2 = (pi +p 2 )/2 and Q = pi - p 2 . 

The data from the L3 collaboration [3] and from the 
CMS colllaboration [5] show that the correlation function 
C'(pi,p 2 ) takes values below unity, contrary to Eq. (3). 
Thus the particles must be correlated and we propose that 
this effect is due to the composite nature of hadrons. 

To implement these space-time correlations, we replace 
formula (1) for the two-particle source function by 

W(p 1 ,p 2 ;xi,x 2 ) = 

w{pi\xi)w{p 2 \x 2 )[l - D{xi - x 2 )\, (5) 

where D{x i — x 2 ) is the cut-off function that satisfies the 
constraint D{x i — x 2 = 0) = 1 and tends to 0 at larger 
distances (above, let us say, 1 fm). Then, the HBT corre¬ 
lation function becomes 
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C (Pi 2 , Q) = C noncoiI (P 12 , Q) - C COII ( PuP2 ), (6) 

today 







where the uncorrelated part C' noncorr (Pi 2 , Q) is given by 
(3), while the correction due to space-time correlations 
reads 

Ceorr = C® + C™. (7) 

where 

„( 0 ) _ J dxidx 2 w(pi-x 1 )w(p 2 ;x 2 )D(xi - x 2 ) 

w( Pl )w( P2 ) ’ 1 j 


C (Q) = 

^corr 

f dxidx2e l( ' Xl ~ X2 ' ,Q 'w(Pi2\x 1 )w(P V 2\x 2 )D(x l - x 2 ) 


One sees that the contribution from the correlation part 
is negative. Moreover, since it obtains contributions from 
a small region of space-time, its dependence on Q is much 
less steep than that of the uncorrelated part. Conse¬ 
quently, at Q large enough C{P\ 2 , Q) may easily fall below 
one. 

To describe the actual measurements one has to take 
into account that particles produced very far from the cen¬ 
ter (e.g. those arising from long-lived resonances) form a 
’’halo” and do not contribute to the HBT correlations [ 8 ]. 
Thus we have 


Coba{Pl1,Q) = 1 -P 2 + P 2 C(P 12 , Q) (10) 


where p 2 is the probability that both particles originate 
from the ’’core”. 

In the ALICE experiment [7] C 0 bs was, in addition, nor¬ 
malized to 1 at some Q o where the influence of quantum 
interference is expected to be negligible. Thus we finally 
have to consider the function 


C 0 bs(Pl2iQ) 


1 -p 2 +p 2 C{P 12 ,Q) 
1-p 2 +p 2 C(P 12 ,Qo)' 


( 11 ) 


Introducing the (measured) intercept parameter A by the 
condition 


l + A = C obs (Pi 2 ,Q = 0) (12) 

one obtains 


C(Pi 2 , Q — 0 ) — C(Pr 2 , Qo ) + A [1 — C(Pi 2 , Qo)] 

(13) 


detail in [ 6 , 9]. In this model, at freeze-out, hadrons are 
created at a fixed (longitudinal) proper time 

T = \Jt 2 — Z 2 = Tf. (14) 

The single-particle source function (in the longitudinal 
c.m.s. system) becomes 

w(p, x) = k 0 cosh 7 ?e - c/cosh,?+v cos< ^ f{r)rdrd(j)dr] (15) 

where ho = Jm 2 + k ^_, whereas 77 , <fi and r are space-time 
rapidity, azimuthal angle and transverse distance from the 
symmetry axis 1 . We have also introduced the notation 

U = /3k 0 cosh 9; V = /3k± sinh#, (16) 

with T = 1//3 being the freeze-out temperature. Finally, 
6 describes the transverse flow by the relation 

sinhd = wr, (17) 


with ui being a parameter. The function f[r) describes the 
transverse profile of the source. 

It was shown in [ 6 ] that the model is flexible enough to 
describe the HBT radii measured by the ALICE collabo¬ 
ration [7]. The function f(r) was taken in the form 

f(r) ~ e-tr-Rf/P (18) 

corresponding to a ’’shell” of the width -J2 6 and radius R. 

Thus the model contains 5 free parameters: T, w, Tf, R 
and <5, which may depend on the multiplicity of the event. 
Their values, giving a good description of the HBT radii 
measured in [7], are given in [ 6 ]. 

4 . Since we treat particles as extended objects produced 
on the hyperbola (14), the longitudinal distance between 
the two hadrons located at the space-time rapidities 771 , 772 
should be calculated along this curve, which yields 

rV2 

d \\ = 

Jr}! 

In the frame where 771 + 772 = 0 we also have t\ = t 2 and 
thus the total distance between particles is 


\/ dz 2 - dt 2 = T f (r ) 2 ~ 771 ). (19) 


d 2 = (x! - x 2 ) 2 + ( 2/1 - 7 / 2 ) 2 + djj = d\ + d 2 . (20) 

Since this expression is invariant under boost in the longi¬ 
tudinal direction, it is also valid in the LCMS system, and 
thus we finally have 

d 2 (x 1 ,x 2 ) =r\ + rl~ 2 nr 2 cos( 0 i - <j> 2 ) + t?(t 7 i - r] 2 ) 2 . 

(21) 


This allows to evaluate the measured correlation function 
in terms of the measured intercept parameter A and the 
evaluated correlation function C{Pi 2l Q). 

Note that in absence of space-time correlations we have 
C{P\ 2 ,Q = 0 ) = 2 , C(P\ 2 , Qo) = 1, and thus p 2 = A, as is 
usually assumed. 

3 . To have an idea on the magnitude of the effect we 
discuss, we have used the blast-wave model described in 


The correlation functions were studied using a gaussian 
cut-off function 

D(x 1 ,x 2 ) = e“ d(xi ’ X2)2/A2 , ( 22 ) 

where A is a constant fixing the scale of the cut-off region. 


1 A11 irrelevant constants are cancelled in the definition of w(p,x). 
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Figure 1: (Color online) Correlation function C Q bs f° r the long di¬ 
rection in the interval 0.2 GeV < Q < 0.8 GeV (normalized to 1 at 
Q= lGeV). The dashed lines describe the results for k± = 163 MeV 
and the two multiplicity classes: N c = 12-16 and N c = 52-151. The 
solid lines describe the results for kj_ = 547 MeV and the same two 
multiplicity classes. 



Figure 2: The same as Fig. 1 but for the side direction. 

5. Using (22), (21) and the source function of the model 
described in Sec. 3, with the parameters taken from [6] 
2 , we have evaluated corrections to the HBT correlation 
functions (6) for all intervals of the particle multiplicity 
and transverse momentum (as measured in [7]), and for all 
three directions of the vector Q. The cut-off distance A « 
2 ry (where ry is the radius of the ’’excluded volume” [10] 
occupied by one pion) was taken to be 1 fm, within the 
range of values given by the earlier analyses [11]. Some of 
the results, obtained using the gaussian D(x i — x 2 ), are 
shown in Figs. 1-3. 

One sees that for the long direction the correlation func¬ 
tion falls below 1 at all multiplicities and transverse mo¬ 
menta of the pair. The depth of the minimum in the long 
direction varies from ~ 0.02 to ~ 0.01 (below 1) when the 


2 The intercept parameter A was taken as A = 0.59 — 0.26 k± 
(where k \ is in GeV) which approximates the data of [7]. 



Figure 3: The same as Fig. 1 but for the out direction. 

HBT radius i?i on g increases from ~ 0.8 to ~ 2 fm. 

In the side and out directions the results are strongly 
dependent on the value of the transverse momentum of 
the pair. At k± < 300 MeV for the side and k± < 400 
MeV for the out direction the correlation functions are 
always larger than 1 in the investigated region. In the side 
direction the correlation function shows a clear structure: 
a minimum followed by a maximum (particularly at low 
multiplicities). At larger k± the minimum below 1 shows 
up in both cases. 

In the side direction the minimum at k± > 300 MeV is 
similar to that found in the long direction. It is about 
twice deeper in the out direction (above 400 MeV). In 
both cases the minimum is deeper when the multiplicity 
increases. Also in this case the change is controlled by the 
corresponding HBT radii. 

To see the sensitivity of these results to the shape of 
the cut-off function D(x\ — x 2 ) we have also considered a 
sharp cut-off which is drastically different from the Gaus¬ 
sian. We have found that the qualitative features are un¬ 
changed, except that the effects of the cut extend to larger 
values of Q. This, however, happens in the region where 
these effects are already small and rather hard to measure. 
Actually, in most cases the results are almost identical 3 , 
provided that the cut-off parameter is ~ 0.75 fm. The only 
exception is the side direction at small multiplicity where 
the difference exceeds slightly 0.02 at Q > 700 MeV. 

We thus conclude that although the shape of the cut-off 
function can influence the details of our results, the general 
qualitative features remain unchanged. 

6. In summary, we have estimated to what extent the 
space-time correlations implied by the excluded volume 
effect modify the HBT correlation functions. 

Our conclusions can be formulated as follows: 

(i) The space-time correlations induced by the finite size 


3 In the relevant region Q > 300 MeV they differ by less than 0.01, 
which is consistent with the expected accuracy of our calculations 
and also with the present experimental accuracy. 
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of hadrons lead to a rich structure of the HBT correlation 
functions, depending on (i) the measurement direction, (ii) 
multiplicity and (iii) the transverse momentum of the pair. 

(ii) The difference between the long and the two other 
directions at small k± is particularly striking. 

(iii) At large k± the minimum below 1 shows up in every 
direction. It is about twice deeper for out than for the long 
and side directions. 

Some comments are in order. 

(i) We have found that the modification of the HBT 
correlation functions are only marginally sensitive to the 
change of shape of the cut-off function D(x i — x 2 ). This 
means that the effect we discuss is, in practice, described 
by a single parameter A. 

(ii) We have been considering the space-time correla¬ 
tions in the source function of two pions, which are a nec¬ 
essary consequence of their composite nature. Naturally, 
there might be also other mechanisms contributing to these 
correlations (e.g. the final state interaction). In this case 
the parameter A should be considered as an effective cut¬ 
off distance which summarizes all contributions. Since our 
calculations show that the measurable effects on the HBT 
correlation functions depend mostly on A (and not on the 
shape of the function D ) it seems hopeless to try to sepa¬ 
rate the various contributions. 

(iii) In our approach the cut-off function is taken in¬ 
dependent of particle density. This approximation seems 
reasonable because, as shown in [6], particle density at 
freeze-out changes only by 10% in the range of multiplic¬ 
ities we consider. Moreover, the dominant effect of the 
changing particle density is expected to be a modification 
of the single-particle source functions of the two pions con¬ 
tributing to interference rather than of their space-time 
correlation described by D{x i — x 2 ). It follows that the 
observable effect of the modification of the cut-off func¬ 
tion due to change of particle density is expected to be 
very small, if any. 

(iv) It is interesting to speculate about the size of the 
effects we discuss in case of heavy ion collisions. Taking 
the source functions in the transverse direction in form of 
Gaussians (which is a reasonable approximation for heavy 
ion collisions) one can easily see that the corrections due 
to finite size of hadrons fall as (A/I?) 2 where R is the 
radius of the system. For PbPb collisions this gives fac¬ 
tor ~ 1/30 compared to the results shown in this paper, 
implying that the expected effects are negligible. Similar 
mechanism should be at work in the longitudinal direction. 
For smaller systems, as those created in p — Pb collisions, 
the effects are also expected to be smaller than in pp. Pre¬ 
cise estimate would, however, require determination of the 
source functions (see [6]). 

Acknowledgements 

We thank Adam Kisiel for useful discussions and for help 
in understanding the ALICE data. Thanks are also due to 
Viktor Begun for an instructive discussion on the excluded 
volume effects. This investigation was supported in part 


by Polish National Science Center grant with decisions No. 
DEC-2013/09/B/ST2/00497 (A.B. and K.Z.) and DEC- 
2012/06/A/ST2/00390 (W.F). 


References 

[1] A. Bialas and K. Zalewski, Phys. Lett. B727 (2013) 182. 

[2] See, e.g., M. A. Lisa, S. Pratt, R. Soltz and U. Wiedemann, 
Ann. Rev. Nucl. Part. Sci. 55, 357 (2005). 

[3] M. Acchiari et al. [L3 coll.] Phys. Lett. B458 (1999) 517; 
P. Achard et al. [L3 coll.], Eur. Phys.J. C71 (2011) 1648. 

[4] G. Abbiendi et al. [OPAL coll.], Eur. Phys. J. C16 (2000) 423; 
P. Abreu et al. [DELPHI coll.] Phys. Lett. B471 (2000) 460; 
A. Heister et al. [ALEPH coll.] Eur. Phys. J. C36 (2004) 460. 

[5] V. Khachatryan et al. [CMS coll.], JEHP 1105 (2011) 029; 

L. Perroni, PoS (WPCF2011) 015. 

[6] A. Bialas, W. Florkowski and K. Zalewski, J. Phys. G: 
Nucl.Part.Phys. 42 (2015) 045001. 

[7] K. Aadmodt et al. (ALICE coll.) Phys. Rev. D84 (2011) 112004. 

[8] P. Grassberger, Nucl. Phys. B120 (1977) 231; M. Gyulassy, 

S. K. Kaufmann and L.W. Wilson, Phys. Rev. C20 (1979) 2276; 
R. Lednicky and M. Podgoretsky, Sov. J. Nucl. Phys. 30 (1979) 
432; T. Csorgo, B. Lorstad and J. Zimanyi, Z.Phys. C71 (1996) 
669. 

[9] A. Bialas, W. Florkowski and K. Zalewski, Acta Phys. Pol. B45 
(2014) 1883. 

[10] R. Hagedorn and J. Rafelski, Phys. Lett. B 97 (1980) 136; 

M. I. Gorenstein, V. K. Petrov and G. M. Zinovev, Phys. Lett. 
B 106 (1981) 327; D. H. Rischke, M. I. Gorenstein, H. Stoecker 
and W. Greiner, Z. Phys. C 51 (1991) 485; G. D. Yen, 
M. I. Gorenstein, W. Greiner and S. N. Yang, Phys. Rev. C 
56 (1997) 2210. 

[11] P. Braun-Munzinger, K. Redlich and J. Stachel, in *R.Hwa (ed.) 
et al. ’’Quark gluon plasma”, p. 491, arXiv:0204013; Y. Hama, 

T. Kodama and O. Socolowski, Jr., Braz. J. Phys. 35 (2005) 
24; L. M. Satarov, M. N. Dmitriev and I. N. Mishustin, Phys. 
Atom. Nucl. 72 (2009) 1390; K. Werner, I. Karpenko, T. Pierog, 
M. Bleicher and K. Mikhailov, Phys. Rev. C 82 (2010) 044904; 
A. V. Merdeev, L. M. Satarov and I. N. Mishustin, Phys. Rev. 
C84 (2011) 014907; A. Andronic, P. Braun-Munzinger and 
J. Stachel, Phys. Lett. B718 (2012) 80; J. Fu, Phys. Lett. B722 
(2013) 144; V. Begun, M. Gazdzicki and M. Gorenstein, Phys. 
Rev. C88 (2013) 024902 and references quoted there. 


4 


